Strong magnetic field in a protoneutron star 



Subrata Pal^, Debades Bandyopadhyay^, and Somenath Chakrabarty^ 
^ Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Calcutta- 700 064, India 
^Department of Physics, University of Kalyani, Kalyani 7^1235, India 
and Inter- University Centre for Astronomy and Astrophysics, Post Bag 4, Ganeshkhind, Pune 411007, India 

We investigate the influence of a strong interior magnetic fleld on the structure and composition 
of a protoneutron star allowing quark-hadron phase transition. In contrast to protoneutron stars 
with noninteracting quark phase, the stars with interacting quark phase have smaller maximum star 
masses than those of neutrino-free stars. The strong field makes the overall equation of state softer 
compared to the field-free case favoring the evolution of a protoneutron star more towards a neutron 
star than a black hole. 
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The understanding of the final journey of a massive 
star, after its fuel has been exhausted is a challenging 
problem . The outcome of it may be a supernova and 
the residue either be a neutron star or a black hole. The 
unique feature about supernova problem is that it in- 
volves all the forces of nature — strong, weak, electro- 
magnetic and gravity. A massive star in its late stage 
of evolution undergoes gravitational core collapse as the 
core exceeds the Chandrasekhar mass. The subsequent 
core bounce occurs when the core density reaches nuclear 
matter value and above. A few milliseconds after the core 
bounce, the hot and lepton rich (neutrino-trapped) core 
settles into hydrostatic equilibrium . The evolution of 
this protoneutron star into a neutron star or a black hole 
completes within a few tens of seconds. 

Observations of pulsars predict large surface magnetic 
fields of B„i 10^^ G The interior magnetic fields are 
a few orders of magnitude larger than the surface fields. 
In fact, the virial theorem [Q predicts large interior field 
of ~ 10^* G or more ||]. One of the plausible expla- 
nations for such a large interior field is that the weak 
field of a progenitor is amplified because of flux conser- 
vation during the gravitational core collapse. The highly 
conducting core results in large ohmic diffusion time so 
that the field is frozen in the core, and consequently not 
manifested at the surface The energy of a charged 

particle changes significantly in the quantum limit if the 
magnetic field is comparable to or above a critical value 
Bm , and the quantum effects are most pronounced when 
the particle moves in the lowest Landau level . In ad- 
dition to this large magnetic field embedded in the dense 
core of the (proto)neutron star, a transition from nuclear 
matter to a stable quark matter is also possible P-pT[|. 

In this letter we investigate the effects of trapped neu- 
trinos on the composition and structure of a protoneu- 
tron star in presence of a strong interior magnetic field 
and also allowing a hadron to quark phase transition in 
the star's interior. The matter inside the protoneutron 
star is highly degenerate and the chemical potential of 
its constituents are a few hundreds of MeV. On the other 



hand, the central temperature of the star is a few tens of 
MeV. Therefore, neglecting finite temperature effect will 
have little infiuence on the gross properties of the star. 
However, the compositional changes caused by trapped 
neutrinos, which are primarily of the electron type Ve, 
may induce relatively larger changes in the maximum 
masses [||. 

We describe the equilibrium conditions of the pure 
quark and baryonic matter and their mixed phase co- 
existing in a uniform background of electrons (e) and 
electron neutrinos (ve) for neutrino-trapped matter in a 
uniform magnetic field Bm along z axis. The neutrino- 
trapped pure quark phase consisting oiu, d and s quarks 
interacting through one-gluon exchange in local charge 
neutral and /9-equilibrium conditions is described by the 
bag model |jl^. The interaction energy density due to 
one-gluon exchange term to order for each flavor with 
Bm 7^ in the zeroth Landau level, 77 = 0, is given by 
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where qf, m*^ and pp are the charge, effective mass and 
Fermi momentum of quark of /th flavor, and p^z — 
Pz[l + Uy /pz]- For the expressions of the Fock contri- 
butions, Uq and Uy , to the single particle energy, we 
refer to Ref. jll|. The QCD coupling constant is defined 
by ac = (?^/47r. The general expression (for all Landau 
levels r/) for the total kinetic energy of the quark phase 
in a magnetic field is 
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where = xO]^^ + \n^{x + 0]^^^ j with 

i = / or e. The notation in Eq. (2) is same as that in Ref. 
^tI , with the first term corresponds to those for quarks 
with df — 3, the second term is that for electrons, and 
the third term corresponds to that for the neutrinos. The 
total energy density of the pure quark phase is then — 

+ E/ '^^/'''"° + Sm + B, where £„, = Bf^/{d,TT) is the 
magnetic field energy density and B is the bag constant. 
The pressure follows fi'om the relation P'' — IJ-fnf + 

fxini — E"^ , where /x/ denotes the quark chemical po- 

(/) / \l/2 

tential, 71/ = (d/g/S„/27r2) X)'^™- g,, _ j 
is the quark density, and I = (e, I'e). The electron den- 
sity is Tie = {eBra/2T:'^)Y^i^o gr|{^J^l-rnl^^Y''^ , and 
He and are the chemical potentials for electrons 
and electron neutrinos. The charge neutrality condition, 
Q'^ = '^f^fnf — He — 0, and the /3-equilibrium con- 
ditions, /id — /i„ -I- /ie — /ii/^ = fJ-s, can be solved self- 
consistently together with the eff'ective masses at a fixed 
baryon number density 71' = (n„ + Ud + ns)/3 to obtain 
the equation of state (EOS) for the deconfined phase. 
For the ease of numerical computation we, however, add 
here the one-gluon exchange term perturbatively to en- 
ergy density and pressure. 

To describe the neutrino-trapped pure hadronic mat- 
ter consisting of neutrons (n), protons (p), electrons (e) 
and Ue, we employ the linear a-u-p model of Ref. 
in the relativistic Hartree approach. The Fock contribu- 
tion to the hadron phase is quite small and there- 
fore neglected. The EOS for this phase is obtained by 
solving self-consistently the effective mass in conjunction 
with the charge neutrality and /3-equilibrium conditions, 
Q'^ = Up— He = and /i„ = fJ'p+IJ'e—fJ'y^ at a fixed baryon 
number density n^. Here Ui and Hi denote the number 
density and chemical potential; the subscript i refers to 
n, p, e and Vf.- The total energy density £^ (given in 
Ref. Inl) and pressure P'' in this phase are related by 

The mixed phase of hadrons and quarks comprising 
of two conserved charges, baryon number and electric 
charge is described following Glendenning The con- 
ditions of global charge neutrality and baryon number 
conservation are imposed through the relations xQ^ + 
(1 - = and = + (1 ^ x)nl, where x 

represents the fractional volume occupied by the hadron 
phase. Furthermore, the mixed phase satisfies the Gibbs' 
phase rules: pp = 2pu + and P^ = P'. The total en- 
ergy density \s £ = x£^ + (1 ~ x)^'^- The neutrino-free 
matter relations can be obtained [ pT[ by putting fj,^^ = 
in the above expressions. 

In the present calculation the values of the dimension- 
less coupling constants for a, lu and p mesons determined 
by reproducing the nuclear matter properties at a satu- 
ration density of uq = 0.16 fm~'^ are adopted from Ref. 
M. The current masses of u and d quarks are taken 



as rriu = m^i = 5 MeV and rus = 150 MeV, and the 
QCD coupling constant is ac = 0.2. Because of trapping 
the numbers of leptons per baryon, Y^e — Ye + , are 
conserved on dynamical time scale. Gravitational col- 
lapse calculations of massive stars indicate that, at the 
onset of trapping, y^e — 0.4. We consider the bag con- 
stant B = 250 MeV fm~"^ which corresponds to the lower 
limit dictated by the requirement that, at low density, 
hadronic matter is the preferred phase. The magnetic 
field, like the baryon density, increases from the surface 
to the center of the star, and consequently the variation 
of the magnetic field Bm with density rib is parametrized 
by the form jl^ 

S,„(nfc/no) = Sr'' + So [1 - exp {-/3(nfc/no)n] , (3) 

where the parameters are chosen to be /3 = 10"''^ and 
7 = 6. The maximum field prevailing at the center 
is taken as Pg = 5 x 10^* G and the surface field is 
P™''^ ~ 10^ G. The number of Landau levels populated 
for a given species is determined by the field Bm and 
baryon density 1^. To the lowest order, the magnetic field 
energy density and magnetic pressure have been treated 
perturbatively. 

A traditional point of view has been that the magnetic 
field is frozen in the neutron star's interior, at least, dur- 
ing the age of the universe, i.e. ^ 10^"^ years, because 
of high electrical conductivity in the core . Only the 
crustal fields, P,„ ~ 10"'^^ G, would diffuse to the sur- 
face and could be observed. However, it was shown by 
Haensel et al. ||T^ that the electrical resistivity of the 
interiors of the star in the normal state (without su- 
perfiuidity and superconductivity) can significantly in- 
crease. This can lead to a dramatic decrease in the inte- 
rior magnetic field decay time which is estimated to be 
t « 2 X 10^^r|P]^'^P|(nb/no)^ years; in the usual nota- 
tion ll^. Soon after the birth of a neutron star when 
the temperature drops below the superconducting criti- 
cal temperature Tc ~ 10^ — 10^*^ (which is model depen- 
dent), it is commonly thought that superfluid protons 
form a type II superconductor in the outer core within 
the density range 0.7no < rib < 2no, and the estimated 
lower and upper critical magnetic fields are respectively 
Pel 10^^ G and Hc2 ^ 3 X 10^^ G (l|]; we have also 
obtained these values in the present model. Of course 
with increasing baryon densities, Ub > 2no, the attrac- 
tive interaction between protons in the ^Sq state (which 
gives rise to proton superconductivity) is diluted, this 
may result in the possible transition to a type I super- 
conductor; the critical field for which is, however, much 
smaller P^ - 10" G 0. 

In our calculation, the magnetic field of 10^* G oc- 
curring at densities « Stiq in the inner core (see Eq. 
(3)) is thus strong enough to suppress superconductivity. 
Therefore, for a neutron star of typical radius P = 10 
km and temperature T = 10^ K, the characteristic diffu- 
sion time for such field, using the above estimate of Ref. 
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I psf , is then f ~ 20 years, and thus should have been ob- 
served at the surface of the neutron star after its birth. 
However, with our choice of the variation of the magnetic 
field with density (see Eq. (3)), at the bulk of the outer 
core (ub ~ (0.7 - 2)no) the field is only lO" - lO^^ G. 
Therefore for T < Tc, the outer core could form a super- 
conducting region. Thus soon after the formation of the 
neutron star, the superconducting outer core region can, 
in principle, trap the much stronger field at the inner core 
and center (in the normal state) preventing its decay to 
the surface. However, the situation becomes much more 
complex in a superconducting state. The homogeneous 
magnetic field then splits into an ensemble of fluxoids 
which contain a strong magnetic field, embedded in 
the field-free superconducting medium for a type II su- 
perconductor. It was also shown that such fields can 
be expelled out of the superconducting region to the sub- 
crustal region in the first ^ lO** years, and subsequently 
they decay to the surface. In any case in what follows, 
for the present study of the evolution of a protoneutron 
star, within the first few tens of seconds when the star 
is hot, into a neutron star or a black hole in presence of 
a strong interior uniform magnetic field, we assume that 
the large fields are frozen in the core and therefore would 
not decay to the surface. 

With the given parameter set of Eq. (3), we show in 
Table I, the mixed phase boundaries for the protoneu- 
tron star (y^e = 0.4) and neutrino-free star {Y^^ = 0) for 
Bjn = and Bq — 5x 10^^ G with interacting quark phase 
(IQP). To examine the measure of the importance of the 
one-gluon exchange contribution, the corresponding re- 
sults for stars with noninteracting quark phase (NQP) 
are given in parentheses. The onset of transition is at 
density ni = uitiq, and the pure quark phase begins at 
n-2 = U2nQ^ where uq = 0.16 fm~'^ is the nuclear mat- 
ter saturation density. In contrast to NQP matter, with 
the inclusion of interaction, the EOS for the quark phase 
becomes softer in the corresponding IQP matter and the 
transition to a mixed phase is delayed. At high density, 
as expected, the EOS for the quark phase is more softer 
resulting in a larger shift in U2 to higher density compared 
to that in ui. As a consequence, the extent of the mixed 
phase is increased in IQP than the corresponding NQP 



matter. For both NQP and IQP matters, neutrino trap- 
ping delays the onset of the transition to higher baryon 
densities and also reduces the extent of the mixed phase 
in contrast to neutrino-free matter for both cases with 
and without magnetic field. 

The maximum masses of the stars Mmax, and their 
central densities ric = UcUq obtained by solving the 
Tolman-Oppenheimer-Volkoff equation are also given in 
Table I for different cases studied. For all cases con- 
sidered, since the central densities Uc are less than U2 
and fall within the mixed phase, the presence of a pure 
quark phase is precluded. With only nucleons and lep- 
tons, without any quark, the maximum masses and cen- 
tral densities of neutrino-trapped(free) stars are respec- 
tively M„iax = 1.699(1. 778)Mq and Uc = 10.400(9.801) 
for Bjn = while those for Bq = 5 x 10^* G are M^ax — 
1.545(1.650)Mo and Uc = 11.075(9.052). Therefore, for 
pure hadronic matter, neutrino trapping generally re- 
duces the maximum mass from that of the neutrino- 
free star, with and also without magnetic field. This 
is caused by the conversion of protons to neutrons, as 
the trapped neutrino leave, which increases the pressure 
more than it decreases by the loss of neutrinos. With the 
introduction of the quarks, which soften the EOS, the 
scenario becomes quite interesting. For the NQP, the 
maximum mass for the neutrino-trapped star is found to 
be larger than that of the neutrino-free star. This re- 
versal in behavior from the pure hadronic case, was in 
fact, first noted by Prakash et al. ||], and an explanation 
sought was that the delayed appearance of the quarks 
in neutrino-trapped star resulted in a stiffer overall EOS. 
However, in contrast to NQP matter, we find for the IQP, 
the trend of the maximum mass is similar to those for 
hadronic stars, i.e., after deleptonization the maximum 
mass of star increases. In order to explain the behavior of 
A/max, we consider the fraction of the gravitational mass 
originating from the pure hadronic part, MhadZ-^max, of 
stars with quark-hadron phase transition. We consider 
first stars in absence of magnetic field. The quarks here 
soften the EOS, thus larger contribution from quarks (i.e. 
small Mhad/-^max) would lead to smaller Mmax- As evi- 
dent from Table I, for stars with NQP, a much larger 



TABLE I. The phase boundaries, ui and U2, and central densities Uc for neutrino-trapped {Ylc ~ 0.4) and neutrino- free 
{Y^^ = 0) stars with maximum masses Mmax/MQ with and without magnetic field that undergo a quark-hadron phase transition 
with interacting quark phase. The fraction of mass originating from the pure hadron phase is Mhad/A^max. The corresponding 
quantities with noninteracting quark phase are shown in parentheses. The variation of magnetic field with density rib is given 
by Eq. (3) with /3 = 10"^, 7 = 6 for Bo = 5 x 10^* G. Calculations are performed for a mean field model of baryons and a bag 
model of quarks with bag constant of i? = 250 MeV fm"''. The nuclear matter saturation density no is 0.16 fm~^. 
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fraction of mass originates from the softer quark-hadron 
mixed phase in neutrino- free stars than neutrino-trapped 
matter. Consequently, M,„ax = 1.610AfQ for neutrino- 
free star with NQP is smaUer than 1.664M0 mass for 
protoneutron star. On the other hand, for stars with IQP 
(with B„i = 0), the pure hadronic part mostly governs 
the EOS (i.e. large Mhad/Afmax) for both neutrino-free 
and trapped matter. As a result the maximum mass here 
follows the trend of pure hadronic stars with B,n = 0, i.e. 
Mniax for Yle = 0.4 is smaller than that for F^^. = 0. This 
also accounts for the fact that the masses of stars with 
IQP are larger than the masses of stars with the corre- 
sponding NQP matter. 

In presence of magnetic field, the hadronic EOS be- 
comes much softer compared to the stiffening caused in 
the quark EOS Therefore, the overall EOS for the 
quark-hadron star with Bm turns out to be softer 
leading to smaller mass than that for the field-free star. 
Thus for Bm 7^ 0, smaller the contribution to the EOS 
from stiff quark part of the mixed phase, smaller would 
be the mass. This, in fact, is revealed by stars with 
NQP for Y^^ = compared to Ylc = 0.4 case. The EOS 
of a star with IQP for Bm ^ is primarily governed 
by the pure hadronic part, their mass is thus reminis- 
cent of that for the pure hadronic stars with _B,„ ^ 0. 
Therefore, the trend of the masses from neutrino-free to 
neutrino-trapped star with magnetic field are same as 
that of field-free stars, both with interacting and nonin- 
teracting quark phase. 

In Fig. 1 we depict a comparison of the composi- 
tion of neutrino-free matter (top panel) and neutrino- 
trapped matter (bottom panel) in a magnetic field of 
Bq = 5 X 10^^ G with IQP. In contrast to Bm = case 
(not shown), the electron fraction, in particular, is en- 
hanced in a magnetic field due to phase space modifi- 
cation H]. In the hadronic sector, the electron fraction 
increases because of neutrino trapping, which in turn in- 
creases the proton fraction due to charge neutrality con- 
dition and thereby reduces the neutron fraction due to 
baryon number conservation. On the other hand, u, d 
and s quarks themselves try to maintain charge neutral- 
ity, resulting in reduction of the electron fraction and 
consequently the neutrino fraction is enhanced in the 
mixed phase. Neutrino trapping also increases the u 
quark abundance in comparison to the neutrino- free star. 

The differences in the abundance of electron, in partic- 
ular within the inner core of different stars investigated, 
would be manifested in the total number of i^e escaping 
during deleptonization. The total number of emitted, 
N^^, from the inner core of mass 0.5Mq during delep- 
tonization from the initial value of Y^e — 0.4 is obtained 
by integrating Ui,^ = n^Y^^ = nh{YLe — Yg) over this re- 
gion of the neutrino-free star. For pure hadronic star, 
N^^ is found to be 9.416 x 10^^ at Bm = 0, while for 
Bq = 5x 10^^ G it is 9.135 x 10^^. On the other hand, 
we obtain for stars with IQP (NQP), the value of N,y^ to 
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Fig. 1. The abundances of neutrino-free (top panel) and 
neutrino-trapped (bottom panel) star with interacting quark- 
hadron phase transition at a magnetic field _Bo = 5 x 10^* G. 



be 1.662 X 10^6(1.994 x lO^^) ^j. ^^^^ ^ q, and for 
Bq — 5 X lO^^G the corresponding values are 1.532 x 
10^6(1.801 X 10^6)^ xhe stars with quark-hadron phase 
transition emit more number of than that by the pure 
hadron stars. This is caused by the drop in the electron 
abundance in the mixed phase as quarks furnish net neg- 
ative charge. The enhancement of electron fraction in 
magnetic field causes a decrease in iV^^ from the cor- 
responding field-free cases. The different scenarios may 
thus be discernible from the difference in the neutrino 
numbers. 

Delayed neutrino emission and possible black hole for- 
mation in the context of SN 1987A have been much de- 
bated issues in recent years. So far there is no observation 
of a pulsar in it. Moreover, the fading away light curve 
leads one to think that it might have collapsed into a low 
mass black hole. Assuming that SN 1987A has gone to 
a black hole, Bethe and Brown ||l^ estimated the gravi- 
tational mass for the compact object in it to be I.56M0 
from Ni production. They argued that any compact ob- 
ject having mass larger than this limit would be a black 
hole. It has been counter argued that the surface mag- 
netic field of a nascent neutron star is weak. It may take a 
few hundred years or more for the magnetic field trapped 
in the crust to reach the surface by ohmic diffusion which 
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would then increase the surface field so that the pulsar 
could be accessible to observation. In the present calcula- 
tion, since the maximum masses of hadron stars increase 
when the trapped neutrinos leave, they will promptly col- 
lapse after core bounce into black holes once the mass of 
1.56M0 is reached, presuming Bethe-Brown limit for the 
maximum neutron star mass. Similar situation arises for 
stars with interacting quark phase (see Table I). In con- 
trast, the maximum masses of stars with noninteracting 
quark phase decrease after deleptonization. Therefore, 
in the event that the maximum masses of these neutrino- 
free stars reach 1.56Af0, the stars will first explode, re- 
turning matter to the galaxy, and then collapse into low 
mass black holes jl9| . We however find in our calculation 
that the maximum star masses both for neutrino-trapped 
and neutrino-free matter in presence of strong magnetic 
field are smaller (than 1.56Mq) than those for the field- 
free stars. Therefore, the presence of a strong interior 
field favors the evolution of a protoneutron star more to- 
wards a neutron star than a black hole even if its mass 
increases after the trapped neutrinos leave. 

Most of the dynamical supernova calculations |p], pO| , pl| 
indicate that a successful prompt supernova explosion 
could be achieved by employing a soft EOS at densities 
larger than nuclear matter density ng. The substantial 
softening of the EOS caused by the strong interior mag- 
netic field (and also by quarks) may provide a viable 
prompt shock mechanism, and thereby merits considera- 
tion in dynamical simulations of explosions. 

In summary, we have investigated the gross properties 
of a protoneutron star in presence of a strong magnetic 
field in the core with a quark-hadron phase transition. 
In contrast to the neutrino-trapped matter with nonin- 
teracting quark phase, the protoneutron star with inter- 
acting quark phase leads to smaller maximum star mass 
compared to neutrino-free case both with and without 
magnetic field. Besides the presence of quarks, the soft- 
ening of the overall EOS caused by magnetic field has 
significant bearing on the evolution of the protoneutron 
star to a neutron star or a low mass black hole. 
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